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Abstract 

Wc obtain the asymptotics for the speed of a particular case of a particle system with 
branching and selection introduced by Berard and Gouere (2010). The proof is based on a 
connection with a supercritical Galton- Watson process censored at a certain level. 

Resume 

Nous etudions un cas particulier de systenie de particules avec branchenient et selection intro- 
duit par Berard et Gouere (2010). Nous obtenons I'asympotique pour la vitesse, en remarquant 
un lien avec un processus de Galton- Watson surcritique censure a un certain niveau. 
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1 Models and Results 

1.1 The censored Galton- Watson process 

For a probability distribution X on non-negative integers, the Galton- Watson process with 

offspring X is the process defined hy Zq = a and 



^k+l — X] ^*^>' 



where the X^^j's are i.i.d. copies of X {Z^_^^ = if = 0). We write (Z/j) = (Z^). We deal here 
with supercritical Galton- Watson processes {i.e. when E(Af) > 1) and to avoid trivialities we impose 
X{0) > 0. In this case, it is well-known (we refer to [Ij for basic facts on branching processes) that 
Zfc dies with a certain probability < q < 1 which is the unique solution in (0, 1) of 

fix) ■.= J2xii)x' = x. 

i>Q 

The first step of this article will be to study a kind of constrained Galton- Watson process, in which 
the constraint is a "roof" at a given height that prevents the process from exploding. 

Definition 1. Given N a positive integer and a probability distribution X on integers, the Galton- 
Watson process censored at level N >2 with offspring X is the process {X^)k>o defined by Xq = N 
and, for k >0, 

mm(^N,ZfIiXk,i) >0, 
otherwise, 



^k+l 



where the Xk^i's are i.i.d. copies of X. 

The process (X^)fc>o is a finite state Markov process which is typically stuck for a long time on 
N, but eventually dies. Let C/tv be the survival time of the Galton- Watson process : 



Un = mm{k > 0,Xf = 0}. 

Let us describe heuristically the asymptotic behavior of Un. When the censored process dies, it 
happens very suddenly: in the uncensored underlying process, the progenies of the N particles pass 
away almost simultaneously in a few generations. This latter event occurs with a probability close 
to , and therefore the censored process is expected to survive a time close to a geometric random 
variable with parameter g^. 

We state this in the following theorem, which may be derived from results of [5]. 

Theorem 1. The following convergence holds in law: 

t/^g^ ^4°° £:(!), (1) 
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where <f (1) stands for the exponential distribution with parameter one. The convergence also holds 
in mean: 



The convergence in law ([T]) is a consequence of (|5],Th.l (3), (4)) and the asymptotic ^ follows 
from ([5|,Th.l (2)). The results of [5] are actually much accurate. As we will only need these simple 
estimates, and for the sake of completeness, we provide in the two following sections an elementary 
and more concise proof of Theorem [T] We also mention |^, in which the Galton- Watson process is 
censored by a function depending on time ; the author obtained a criterion for the degeneracy of 
the process. 

1.2 A connection with a branching-selection process 

The present authors considered Theorem [T] when trying to find an asymptotic of the speed of 
a very particular case of a certain branching particle system with selection, studied by Berard and 
Gouere ([2], Section 7, Theorem 5). The branching-selection process we will study is a generalization 
of their Bernoulli branching-selection process and can be described as follows. For a distribution 
on such that the expectation of X is strictly greater than 1 and X + X' > 1 a.s., the 
particle system is the discrete-time particle system of N particles moving on Z and starting from 
the origin such that, at each time unit, 

1. Each of the particles is replaced by X particles just on the right of that particle and by X' 
particles on the same position, independently of each other ; 

2. Among all these new particles, we keep only the N rightmost particles. 

It is convenient to see the process of the locations of the A^ particles as a finite point measure. 
We set Yf(^ = N5o and write 



where Yi^{i) is the number of particles at time k at the position £ ; by construction we have at any 




(2) 



time Ef>o^/f (^) = ^- We write 



maxY"^ := max 



{^>0;y,^(^)>0} 
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for the position of the rightmost particle at time k. We have maxY"^ < by construction. Since 
E(Af) > 1, and since Yj^-^^ > by the property of Af' , it is Ukely for N large that maxY^^ is close 
to k. Remark also that in this process, all the particles at time k are on {maxY"^^ — l,maxY^^} 
after the selection. 

Berard and Gouere in fact studied the particle system defined as follows. At every time unit each 
particle is duplicated and moves one step forward with probability a G (5,1) and stays put with 
probability 1 — a. With our notations, this corresponds to the case where X is a Bernoulli variable 
B{2,a) with a > 5, and X' = 2-X (we will denote their model as the Bernoulli branching-selection 
process). They proved that (maxY^^)/A; converges almost surely to a constant f7v(a), and noticed 
that vn{ol) lies between (1 — exp(— ci(a)A^)) and (1 — exp(— C2(a)A^)) for some positive constants 
ci(a) and 02(0). The stress of their paper was a large class of distributions for the walk performed 
by the particles, and the case of the Bernoulli random walk with a > 1/2 was in fact just mentioned 
as a degenerate case. The existence of the asymptotic speed VN{a) would apply as well in our case, 
but not their proof for the bounds on this speed. 

The key point for our result is that the process of the number of particles that are at the rightmost 
possible position {Yj^{k))k has the same law as the censored Galton- Watson process {X^)f^, when 
the offspring distribution is X. 



Theorem 2. For the branching-selection process with the distribution (X,X') on such that 
K(X) > 1 and X -\- X' > 1 a.s., we have 



lim 

Af-5-oo 



1 - VN{a) 



= 1 



where q is the extinction probability of a Galton-Watson process whose offspring is X. 



In particular, for the Bernoulli branching-selection model studied in [2j with a > 



1 

2' 



lim 

Af-s>oo 



1 - VN{a) 



= 1 



(.Qa) 



where 



Qa = 



1 - 2a(l - a) - y^l -4a(l - a) 



2a2 
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is the extinction probability of a Galton- Watson process whose offspring is the Binomial distribution 
with parameters (2, a). 

Remark that the exact law of X' has no influence, and we could have always taken 

^' = ix=o- 

Using A" allows us essentially to have a true generalization of the Bernoulli branching-selection 
process. 

The article is organized as follows. The proof of Theorem [T] is decomposed over Sections [2] and 
|3] A key ingredient in Section [3] is to compare our censored process at the last time it reaches the 
level N, and the classical Galton- Watson process starting at N and conditioned to die. The main 
contribution is the application to the branching-selection particle system and is given in Section |4] 

2 Preliminaries 

We will repeatedly use the following left-tail bound for a sum of copies of X : 

Lemma 1. Let Xi, . . . ,Xjv be i.i.d. copies of law X. There exist b,c > such that, when N is 
large enough, 

P{Xi + --- + Xn <N{1 + b)) < exp(-cA^). 

To prove this, we choose M such that E[min{Xi, M}] > 1. We then can apply ([4J, III Chap. 4) 
to the bounded random variables min{Xj,M} and we get the desired bound. This implies in 
particular that 

P(Xf+i <N\X^ = N)< eM-cN). (3) 

Note that throughout the paper, c stands for a generic positive constant, which might differ at 
each appearance. We introduce the last time for which our process is equal to A^. This variable will 
turn out to be equivalent to Un and more simple to approximate by a geometric variable. 

Definition 2. Let 

Vn = max{k > 0, = A^}, 
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and let be the probability that {X^)f^ does not ever hit N after time zero: 

QN = <N,yk>0)= F{Vn = 0). 

We deduce from Q that 

QN < P(Xf <N)< exp{-cN). (4) 

We will compare qjsf with the probability that the classical Galton- Watson process (Z^) start- 
ing from N particles dies out. By independence of the progeny of the particles, we have 
P((Z^) dies out ) = . By conditioning over the first passage over some level n, we get 

P((Zfc) dies out | 3A; > such that Zk > n) < g". (5) 



3 Proof of theorem [T] 

The main tool of the proof is the following lemma, which says that the Galton- Watson process 
(Z^), conditioned to die, roughly behaves like after its last passage in A^. Recall that q stands 
for the probability that Z dies while q]\i is the probability that does not ever hit N. 

Lemma 2. There exists c > such that for N large enough 



qN 



< exp(-cA^). 



Proof By Lemma [T] we can take 6 > and c > such that 



P(Zf < (1 + b)N) < exp{-cN). 



The point is that if {Z^)f^ dies out, then with high probability it does not hit after time fc = 0. 
Let R = inf {A; > 1, Z^ > iV} (with the convention that R = +oo is does not hit N after time 
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zero) , 



F{R < oo; dies out) 

< P(Zf < (1 + b)N; R<oo; dies out) + P(Zf > (1 + b)N; Z^ dies out) 

< P(Zf < (1 + b)N; R<oo; Z^ dies out) + g(i+^)^ 

< Yl ^(^R R< +oo; Z^ < (1 + b)N; Z^ dies out) + q^^+^^^ 
e>N 

< J2 ^(-^^ d.o.| Z^ = e;R< +oo; < (1 + 6)Ar)P(Zij = i-R< +oo; < (1 + 6)A^) 



1>N 

Jl+h)N 



< J2 ^(^^ d.o.| = ^; < +oo)P(Zf < (1 + b)N) + q^^+^^^ 



e>N 

i>N ^ 

Finally, for N large enough 

P(Z^ dies out) - F{Vn = 0) = P(Z^ dies out, R < oo) 

< g^exp(-c'Ar), (6) 

which finishes the proof. □ 

By construction Un > Vn, we first prove that Un is close to Vn, in the following sense. 
Lemma 3. The sequence (iq?^)iv>i converges to 1 in probability. 
Proof. By definition we have Un > Vn + For e > 0, 

P (y+V^ -!<-£)< HUn - {Vn + 1) > VeN) + ¥(uN < . 
We study the first term using a standard technique in branching processes theory: 

F{{Zl) dies;ZK>0) = q-fK{0), 
where fx is both the K-th iterate of / and the generating function of Z^. Our assumptions imply 



that \q - /k(0)| < exp(-cK) for some c > ([T],Th.l,chap.I.ll). 

By the Markov property, the process {X^_^^)i>o conditioned to = N has the same law as 
(X^)£>0, it follows that 

- (Vn + 1)>K) = F{Un > K + 1\ Vn = 0). 

Recall (|6j) : 

F{Vn = 0) > P(Z^ dies out) - 2q^ exp(-c iV) > /2 



for N large enough. We then write 

F{Un > K + 1\Vn = 0) = 

< 

< 
< 

where we finally used Sec. I.ll). We obtain 

F{Un -Vn-1> VeN) < CiVexp(-cVeiV) 



P(Zj^+i >0\Z^ <N for each 
P(Z^+i > I (Zf ) dies out ) X 

NF{Zl-^^ > I (Zf ) dies out ) 
2^exp(-cK) 



A: > 1) 

P ((Zf ) dies out ) 
P(Zf < N for each A; > 1) 



(7) 



for some C and c > and for any A'^ > 1 and e > 0. 

The second term is handled thanks to Q as follows. For N > 1/e, 



f/«<L^ivj) < p(v«<L^iVj 

< P < [^iVj, V/V = V 



< f(^3v< = N and X^^^ < 

< ^A^exp(-cAr) 
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for a certain c > 0. Letting N go to infinity finishes the proof. □ 

We introduce another process and an associated time whose properly renormahzed law will 
converge and which is equivalent to V/v- Let A/^ be the fc-th passage in N of the censored process: 
Aq = and for A; > 0, A^+i = inf{£ > A^, = N}, with the usual convention that inf(0) = +oo. 
Let T be the survival time of A : 

T = sup{k > 0, < oo}. 

It is clear by contruction that T has the law of G^Qn) — 1, where G{r) is a geometric of parameter 
r. By the Markov property of our branching process, the variables Ak+i — A^ are independent and 
identically distributed. With high probability, A^+i = Ak + 1, we can be more precise: 

Lemma 4. There exists c > such that for all N 

¥{Ak+i >Ak + l)< exp(-ciV). 

There exists c,C > such that for all K > 1, for all N large enough, 

¥{Ak + K < Ak+i < oo) < Cexp ("C^^ - cN^ ■ 

Proof. The first assertion is just ([3|. The second one is a consequence of the following inequality: 
for any < i < N 

P(0 <Zl<N for all k < K) < C exp ("Cj^^^ • (8) 

To prove so, take 1 < /i^. < /i < +oo. The sequence Zf^//j,^ tends to infinity with positive probability 
(see ^j,Th.3 chap.LlO). Hence we have positive constants 6,13 such that, for any integers k,n, 



> Si,':) > ml > Sf,':) > p. (9) 

Now, set m = [log(A^/(5)/ log(^^)] . We assume that N is large enough, such that m>l. We first 
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assume that K > m. 



{0<Zk<N for all k < K} C {Z„, < iV} n {Z^ > and < A^} 

n {Z2m > and < A^} Pi • • • n {Zi^K/m\m-m > and Z^K/mlm < N} . 

Thanks to ([9|, the probability of the right-hand side is smaller than (1 — li K < m then 

K/logN < 2/log/i* for N large enough. By choosing C large enough, the right-hand side in ([8| is 
greater than one, and the claimed inequality also holds. 

To finish the proof of the lemma, it suffices to remark that for K > 1, 

F{Ak + K < Ak+i < oo) < ¥{Ak + 1 < A^+i < oo)P(0 < < iVV^ < K) 

< C exp(—cN) exp —c 



logiV . 

□ 

Lemma 5. The sequence converges to one in probability. 

Proof. By construction we have T < V/y. In order to bound the difference, we first notice that, T 
being a Q^Qn) — 1, we have for large enough the bound 

P(r < 1/^) < (1/^ + l)P(r = 0) < (l/V^ + l)qN < 2^ < 3q^/^ (10) 

where we finally used Lemma [2] We now write 

k=0 k=0 

where the (Ak) is a sequence of i.i.d. random variables with the same law as yli — 1 conditioned by 
^1 < oo, and the sequence is independent from T. 
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Now take e > 0. By Lemma [i] for N large enough we have E(^i) < e/2. 
P (y,v + 1 - (T + 1) > e (T + 1)) = P ("^ Ifc > e{T + 1)^ 

\fc=0 / 

< P ("^ Ifc - E{Ak) > (T + l)e/2\ 

\k=0 / 

< P ("^ Ifc - E(Ifc) > e{T + l)/2; T > 1/^^] + P (T < 1/^) 

\fc=o / 

^ Yl IP fe^fc-IE(^fc) >e(i + l)/2) +IP(r< l/V?^) 
^ Yl pfe^-E(Ifc) >e(t + l)/2j +P(T< l/Vg]^), 

where we finally used Lemma |2] The second term on the right-hand side goes to zero thanks to 
(10). Let us now handle the sum. Using Lemma |4] we have 

P (^fc >K^= ^{K <Ai< oo)P(Ai < oo)"^ < 2exp(-cK/logiV) 

for a certain c > 0. From such a tail, we deduce (see sec.IIL4) that there exists c > such that 
for all t, e and N, 



/ 



and the sum goes to zero. □ 
End of proof of Theorem [7J We now write 

where, when N goes to infinity, 

C/Ar/(Viv + l) —^1 (in prob.) by Lemma jsj 

/qn —^1 by Lemma 

qn / {T + 1) — )• <f (1) (in law) since T is geometric with parameter qi\f, 



{Vn + l)/(T +1) -^1 (in prob.) by Lemma [5] 
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This proves the convergence in distribution. For the convergence in mean, recall that one may write 

T 

¥.[Vn] = E[r + ^ Ifc] = E[r] (l + E[Ii]) = + o(l)), 

fc=0 

thanks to Lemma |2] and |4] Since E[C/Af] = IE[V/v] + '^\Un — ^7v]) we finally obtain with ([T]) 

nuN] ~ nvN] ~ {i/qf. (11) 

which will be useful in the next section. □ 

4 Application: Branching-selection particle system 

We now describe more precisely the connection with the branching-selection process defined in 
the introduction, in order to exploit the results of the previous section. Recall that {(.) is the 
number of particles at time k at the position I and that vjsi := limfc_j.oo(maxyj|^)/A; (a.s.). 

The connection with the censored Galton- Watson Process in is the following Lemma, whose 
proof is immediate by construction of the process. 

Lemma 6. In the branching- selection process with a pair {X^ X') of laws on the number of 
particles that are at the rightmost possible position has the law of a censored Galton- Watson process, 
when the offspring distribution is X : 

In particular, if we define as the last time at which the N particles are at the rightmost 
possible location: 

= ma^ {k;Y^{k) = N} , 

then ^= V/v, where V^, defined in the previous section, is the last time at which the censored 
Galton- Watson hits N. 

We now are ready to estimate v^. 
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Figure 1: A realization of the true process Y with = 3 particles, for the particular case of the 
Bernoulli branching-selection process. At each step are represented the 6 particles (after replication) 
and the 3 leftmost particles are striked. Here, we see that Vi = 3, since from time k = 5 there can 
be no more particles at the rightmost possible position, and because at time 4 only one particle is 
on the rightmost possible position. 



Lemma 7. For all integer N > 1, 

where the progeny distribution defining Vn is X . 

Proof. The proof is inspired by that of Proposition 4 in f2] . The main idea is to design a dominated 
process Yk moving slower than the true process Y^. 

Let us skip the exponent N in order to lighten the notations, and consider an i.i.d. family 
{^i.k,ii'^iki)t>o.k>Q,i&[i,N] of integer variables with the same law as {X,X'). 

We sample the particle system (1^) by means of these variables: the population 1^ at time k 
being given, we let 

00 ii-W 

Tk+i = {Xi^k,ik+i + ^lk,i^i) 

1=0 i=l 

be the locations of the particles after branching. The population l^+i is then obtained by keeping 

only the N rightmost particles among Tk-\-i- 

We modify Y in the following way. At time 1/^ + 1, let continue the process as if all the particles 

were at position at time + 1. More precisely, set (Y^^^) = (Yk) and introduce a process 

(yP^) starting from 
V /fc>yi 

Y^';>^^ = N6yr. 
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Figure 2: A realization of the associated dominated process: at time + 1 = 4, we let the process 
restart with all particles at position 3. 



Then, ior k > + 1, Y^'^^ is sampled from Y^^^ exactly as Ffc+i is sampled from Y^ (with the 
same (A^^fc^j, A"^' ^ ,.)'s). The main point is that, at time + 1, the N particles of Y^^^-^ are at 
positions or + 1 while the particles of ^yi^^^^ are at position V^. Hence, the point measure 
^yi+i dominates Yy^^^, and this domination will continue throughout the process, since the same 
{X,X'ys are used to generate Y'^^^ and Y^'^\ 

Now, let V'^ be such that + 1 + V'^ is the last time k at which the N particles are at the 
rightmost possible location (which is position k — 1) for y^^^: 



= max jfe; Y^'^\k - 1) = A^} - {V^ + 1) 



The random time V"^ has the same law as and, as a function of ^Xj^^k,i,'^eki'i^ > ^^jj is 
independent of V^. We define recursively similar processes Y^^\Y^^\ . . . and random variables 
V^,V'^,.... Setting also 

po ^ pi+i ^Y' + V' + l, 

the r*'s are renewals. 

Let us introduce a last process {Yi^)k>o as follows: if fc G [r*~^,r*), then 



n = Y^^ 



One can see in Fig. 1 the first steps of a realization of the process (1^), and in Fig. 2 the corresponding 
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modificated process (Ifc). By construction, Yj, dominates for each k, and in particular maxYj. > 
maxYk- Let us also note that at each renewal F* the quantity maxl^ is decreased by one. Set 
Ik be the renewal process associated to the renewals F*, that is 1^ is the only integer i such that 
r* < < r'+^. On the one hand we have 



-pi+l pi ^1 _j_ 



SO by applying the renewal theorem (see for instance [3] Chap. 3. 4) to the renewals r^,r^, . . . , we 
get 



On the second hand we have 



lim -E(4) ^ -^-4, . 



max Yk = k — Ik- 



Hence we get 



^-E[maxyfe] > ^E[maxyfc] 1 ^ 



k ^ ^ - k ^ ^ E[V^] + 1' 

To conclude, recall that thanks to the connection with the censored Galton- Watson process 
E[yi] =E[V/v]. □ 

Lemma 8. For all integer N > 1, 

VN <1 



E{Un) 

Proof. The proof is an adaptation of the previous one, replacing by C/^, which is the first time 
at which there is no more particle at the rightmost possible position: 

= mm{k;Yi^{k) = 0} , 

(On the example drawn in Fig. 1, one has = 5.) Then, has the same law as the defined 

in Section [2] In a similar manner to the previous proof, we let the process restart at time as if 

all the particles were at position — 1. The end of the proof is similar. □ 

We now combine these two estimates oi with the results of the previous section in order to 
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prove Theorem |2] : 



1 1 

<1-VN < 



E{Un) - E{Vn) + 1 

and both sides, thanks to (11), are equivalent to . 
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